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In view of the recent WA98 data of 7r° spectra from central Pb + Pb collisions at the 

CERN SPS, we analyze the production of neutral pions for A + A collisions across the 
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periodic table at yfs = 17 AGeV and 200 AGeV within the framework of the parton- 
cascade model for relativistic heavy ion collisions. The multiplicity of the pions (having 
Pt > 0.5 GeV) in the central rapidity region, is seen to scale as ~ (N part ) a , where N part 
is the number of participating nucleons, which we have approximated as 2A for central 
collisions of identical nuclei. We argue that the deviation of a (~ 1.2) from unity may 
have its origin in the multiple scattering suffered by the partons. We also find that the 
constant of proportionality in the above scaling relation increases substantially in going 
from SPS to RHIC energies. This would imply that the (semi)hard partonic activity 
becomes a much cleaner signal above the soft particle production at the higher energy of 
RHIC, and thus much less dependent on the (lack of) understanding of the underlying 
soft physics background. 



Recently, the WA98 Collaboration has published |IJ data for the production of neutral 
pions up to transverse momenta of p± ~ 4 GeV/c, in central Pb + Pb collisions at 160 
A GeV/c incident momentum, corresponding to \/s ~ 17 A GeV. Two most interesting 
features of these data emerge when compared to corresponding data from pp collisions 
and collisions involving lighter nuclei [[J: a) an approximate invariance of the spectral 
shapes, i.e., a near indepence of the slope of the neutral pion p± spectra; b) a simple 
scaling of the 7r° with the number of participating nucleons, if the number of participants 
is large (£30). In the present work, we use the event generator VNI which embodies 
the physics of the parton-cascade model Q for ultra-relativistic heavy-ion collisions to 
analyze these observations. The model attempts to describe the nuclear dynamics on 
the microscopic level of particle transport and interactions, by evolving the multi-particle 
system in space-time from the instant of nuclear overlap all the way to the final-state 
hadron yield. For details we refer the interested reader to Refs. 0|§. 

A simple consideration may illustrate the features of particle production within this 



approach and its relevance to multiple parton scattering. Let x denote the number of 
partons in each nucleon, and let each parton suffer v collisions during the partonic stage. 
Assuming that each virtual parton radiates r partons, we see that the number of produced 
partons will vary as 

^"partons OC U (I + r) X A , (1) 

and as we expect, -Whadrons oc -/Vp ar t ns , one realizes immediately that if the partons 
interact only once, the multiplicity of the partons, and hence the multiplicity of hadrons, 
will scale as A. It is also clear that if every parton interacts with with every other parton 
then v oc A, and the number of materialized partons would scale as A 2 . That can happen, 
if the system would live for an infinitely long time. However, this is not the case. Rather 
than that, in relativistic heavy ion collisions, the partonic matter will expand, dilute, and 
eventually convert into hadrons. Thus a given parton may undergo v ~ R/X interactions; 
where R is the transverse size of the system and A is the mean free path of the parton. 
Noting that R ~ A 1 ^ 3 , we immediately see that the number of materialized partons, 
and hence the number of produced particles would scale as ~ A^ 3 . An experimental 
verification of this scaling behaviour could be a direct manifestation the formation of a 
dense partonic matter! 

We shall demonstrate now that these simple considerations are indeed confirmed by a 
detailed simulation with the event generator VNI on the basis of the parton-cascade/cluster- 
hadronization model. We first consider the recently measured transverse momentum dis- 
tribution of 7r°-production in central collsions of Pb + Pb at CERN SPS obtained by the 
WA98 collaboration [I]. To make contact with the experimental data, the simulations 
were done for the range of impact parameters < b < 4.5 fm, which corresponds to 
10% of minimum-bias cross-section. The result of our model calculation, shown as the 
solid histogram in Fig. la, is seen to be in decent agreement with the experimental mea- 
surements. The model results do not include the final-state interaction among produced 
hadrons yet, but it is likely 0] that the agreement will further improve once the effect of 
cascading hadrons is included. The dashed histogram in Fig. 1 gives the soft contribu- 
tion while the solid curve gives a hydrodynamic prediction (without the contribution of 
resonance decays). 

In Fig. 2 we plot our results for the p± spectra of 7r°'s for a number of central A A 
collisions at y/s = 17 A GeV for various A + A systems from A = 16 to A = 197. One 
observes that they are almost identical in shape with a universal slope for p± <, 1.5 GeV/c. 
On the other hand, the deviations appearing at larger p± for heavier systems are indicative 
of enhanced multiple scattering there. Similar results (not displayed here) were obtained 
at RHIC energies. In order to verify this scaling more closely, we have calculated, as a 
function of the nuclear mass number A, the production of 7r°'s in the central rapidity 
region (—0.5 < y < 0.5) having transverse momenta p± > 0.5 GeV/c. The latter choice 
minimizes the influence of pions having their origin in decay of resonances. This kinematic 
window was motivated [|I[] by the WA98 collaboration in their measurement of the 7r° yield. 
Fig. 3 displays the simulation results for central A + A collisions across the periodic table, 
at CERN SPS center-of-mass energy y/s = 17 A GeV, while Fig. 4 shows the same for 
RHIC energy yfs = 200 A GeV. The solid lines are fits to the model results, represented 
by the symbols, and scale as 
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Figure 1. Transverse mass spec- 
tra of neutral pions in central col- 
lision of 158 AGeV Pb + Pb colli- 
sions. The solid histogram repre- 
sents our result from the parton- 
cascade/cluster-hadronization 
model. The dashed histogram 
gives the contribution of the soft- 
part. 
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Figure 2. Transverse mo- 

mentum spectra of neutral pi- 
ons in central collisions of iden- 
tical nuclei at E cm = 17 A GeV. 
The symbols are our results for 
various systems, from O + O up 
to Au + Au. All results are nor- 
malized to the case A = 16. 



N^o ex (N part ) a , (2) 

where N part = 2A is the number of participating nucleons, and a being extracted as: 

(3) 



a 



1.16 at y/s = 17 A GeV 
1.23 at v/i = 200A GeV ' 



1.2 is in excellent agreement with the corrected WA98 



It is interesting to note that a 
results Q. 

For comparison, the dashed lines correspond to a linear scaling ~ {Npart) 1 ' , whereas 
the dashed-dotted lines indicate a hypothetical scaling with the nuclear overlap factor 
Taa(P = 0) ~ (Npart) 1 A2 . A linear scaling would reflect a single-collision situation, and 
a scaling with Taa would indicate a Glauber-type multiple-collision scenario, in which 
nucleons suffer several collisions along their incident straight-line trajectory, without de- 
flection and without energy loss. Comparing the three curves, we can conclude that our 
simulation results rise significantly slower than with Taa, because firstly, the particles 
change direction through the collisions, and secondly, they are subject to a collision time 
of the order of the inverse momentum transfer, during which they cannot rescatter. On 
the other hand, our calculated tc° yields grow much faster than linear with A, due to 
multiple scatterings. 
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Figure 3. Mass number scal- 
ing of pion production in the 
central rapidity region at CERN 
SPS energies. The symbols rep- 
resent the results of our simula- 
tions, and the solid curve is a fit 
to these results. 



Figure 4. Same as above for 
BNL RHIC energies. Both the 
increased constant of proportion- 
ality as well the exponent for 
scaling imply an increased mul- 
tiple scattering at the (higher) 
RHIC energies. 



In summary, we have demonstrated here that the observed scaling of the number of 
produced particles with the number of participants, in heavy-ion A + A collisions, as well 
as the approximate shape-independence of the transverse momentum spectra, are satis- 
factorily reproduced by the parton-cascade / cluster-hadronization model. We must add 
that a more accurate description of the p? spectra at larger p? will need a readjustment 
of the hadronization model. 

Dr. Klaus Geiger, a dear colleague, a very good friend, and a brilliant physicist was 
killed in an air-crash in September 1998. This is one of the last pieces of work which we 
completed together /[J/. 
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